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In this paper we analyze the asymptotic behavior of Cardassian cosmological models filled with 
a perfect fluid and a scalar field with an exponential potential. Cardassian cosmologies arise from 
modifications of the Friedmann equation, and among the different proposals within that framework 
we will choose those of the form 3 H 2 — p oc p n with n < 1. We construct two three dimensional 
dynamical systems arising from the evolution equations, respectively adapted for studying the high 
and low energy limits. Using standard dynamical systems techniques we find the fixed points and 
characterize the solutions they represent. We pay especial attention to the properties inherent to 
the modifications and compare with the (standard) unmodified scenario. Among other interesting 
results, we find there are no late-time tracking attractors. 

PACS numbers: 98.80.-k, 98.80. Jk. 


I. INTRODUCTION 

Observations indicate that type la high redshift super¬ 
novae (SNIa) are dimmer than expected (see for most 
recent results), and the mainstream interpretation of this 
result is that the universe is currently undergoing accel¬ 
erated expansion. Further observations, like those of the 
Cosmic Microwave Background (CMB) (2j or Large Scale 
structures (LSS) §j, suggest that two thirds of the energy 
density of the universe correspond to dark energy. 

Even though the most popular explanation to the late- 
tirne acceleration in the universe is the existence of some 
kind of dark energy (perhaps a scalar field), this is not the 
only possibility. Recently, Freese and Lewis pj proposed 
the so-called Cardassian models as an alternative expla¬ 
nation which involves only matter and radiation and does 
not invoke either vacuum energy or a cosmological con¬ 
stant. In these models the universe has a flat geometry, 
as required by measurements of the cosmic background 
radiation [j|, and it is filled only with radiation and mat¬ 
ter (baryonic or not). The Friedman equation is modified 
with respect to its usual form by the addition of a term 
in its right hand side, specifically 

3 H 2 = p + ap n , (1.1) 

in units such that 87r/m^ = 1 and with a > 0 being an 
arbitrary constant. 

For n < 1 the second term becomes important if 
z < 0(1). From there on it dominates the Friedmann 
equation and yields a oc t 2 / 3rl for ordinary matter, so 
there will be acceleration provided n < 2/3. There 
are two main (possibly unrelated) motivations for the 
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p n modifications: (1) As shown in Jj|, terms of that 
form typically appear in the Friedman equation when 
the universe is embedded as a three-dimensional surface 
(3-brane) in higher dimensions. (2) Alternatively, these 
functions may appear in a purely four dimensional the¬ 
ory in which the modified right hand side of the Friedman 
equation is due to an extra contribution to the total en¬ 
ergy density. One will then regard the right hand side of 
the Friedman equation as corresponding to a single fluid, 
with an extra contribution to the energy-density tensor 
in the (ordinary four dimensional) Einstein equations. 

The interpretation of the Cardassian expansion as due 
to an interacting dark matter fluid with negative pressure 
was developed in [?|. The Cardassian term on the right 
hand side of the Friedman equation is interpreted as the 
interacting term and gives rise to the effective negative 
pressure which drives the cosmological acceleration. 

Interestingly, Cardassian models survive several obser¬ 
vational tests, the most significant being that it allows for 
a universe consisting of just matter and radiation. The 
energy density giving a closed universe p c is much smaller 
that its counterpart standard cosmology p qo \d (specifi¬ 
cally p c = ^ old [l + (1 + z e q) 3 ( 1 - n)] 1 ), and matter 
alone is enough to provide a flat geometry. 

With regard to the observational tests that depend 
only on the scale or the Hubble factor, in the late¬ 
time regime Cardassian models filled with just matter 
(p oc a~ 3 ) are indistinguishable from perfect fluid models 
with a p = (7 — 1 )p equation of state under the identifica¬ 
tion n = 7. These perfect fluid models are in turn kine¬ 
matically equivalent to scalar field (quintessence) models 
with an exponential potential. In this way, the Cardas¬ 
sian model can make contact with quintessence with re¬ 
gard to observational tests. Interestingly, observational 
tests seem to favor n < 0, so that as ymp totically one 
would get a phantom equation of state |19j . 

However, the equivalence between Cardassian and per- 
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feet fluid models is not extensible to the dynamical realm, 
the evolution of perturbations may differ significantly, 
and this can lead to discrepancies for instance in ob¬ 
servational tests associated with the cosmic microwave 
bakground. 

Nevertheless, as stated in Jj, questions of interpreta¬ 
tion remain open, because in the Cardassian model mat¬ 
ter alone is responsible for the accelerated behavior, and 
yet the universe can be flat. 

Here we perform a dynamical systems analysis of Car¬ 
dassian models. These techniques have been used for 
exploring Cardassian models filled with baryonic matter 
only j 22 j, but we give one step further by allowing as well 
for a scalar field component (non-baryonic matter). We 
investigate for the early and late time regime of the mod¬ 
els the existence of tracking solutions, i.e., we look for 
solutions representing universes with non-negligible and 
proportional fractions of both baryonic and non-baryonic 
matter. Tracking solutions are particularly interesting 
because their dynamical effects mimic a decaying cos¬ 
mological constant (see 00 for seminal references). 
Such solutions would be devoid of the fine-tuning prob¬ 
lems posed by a cosmological constant precisely because 
of the independence on the initial conditions. 

Using standard dynamical systems techniques we will 
investigate the interplay between the modifications to the 
Friedmann equations and the features and existence con¬ 
ditions of the mentioned tracking solutions. In this first 
approach to the problem we will choose an exponential 
self-interaction potential 000 for the scalar field 
because only those potentials allow using the evolution 
equations for constructing an autonomous system of as 
low a dimension as possible. As shown in [Uj; for an y 
other potential the number of dimensions will be higher 
if the system is to remain autonomous. In that refer¬ 
ence the existence of tracking Cardassian cosmologies was 
shown within a more general context, but their analysis 
is not as detailed as ours. 

The remark has to be done that even though obser¬ 
vations seem to favor n < 0 J 2 l|, we will restrict our 
numerical examples to cases with 0<n<2/3so that we 
can establish comparisons between our results and those 
in 0 - Nevertheless, most aspects of our analysis do not 
depend on the sign of n, but just on the requirement that 
n < 1 . 

The plot of the paper is as follows. In Section II we 
study the phase-space, find its fixed points and charac¬ 
terize them. In Section III we discuss the cosmological 
consequences of the attractor solutions. Finally, in Sec¬ 
tion IV we outline our main conclusions, including our 
main result that there are no tracking late-time attrac¬ 
tors. 


II. PHASE-SPACE 

The evolution equations for a flat Friedmann- 
Robertson-Walker (FRW) Cardassian cosmological 


model filled with a scalar field <fi with self-interaction 
potential V(<f>) = £ exp(— %/ 6 A^>) and a barotropic perfect 
fluid with equation of state p 1 = (7 — 1 )p 7 are 


2 H 
0 + 


[iPi + 4> 2 ^j 


'YP'y+4’ 

dV{V) 
d<j) 


(1 +nap? ot 1 ) = 0 
= 0 , 


p 7 + 3y Hpj = 0, 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


where for the total energy density ptot we have 

Ptot = 2^ 2 + ^(0) + p 7 - (2.4) 

The evolution equations are in turn subject 

to the constraint 


H 2 = ip tot (1 + ap?~V) . (2.5) 


Here and throughout cr, k and A will be free parameters, 
and we will restrict ourselves to the n < 1 case. 

Experience has demonstrated that dynamical systems 
methods can be used to describe the evolution of cosmo¬ 
logical models by means of past and future attractors. In 
order to cast our set of equations as a dynamical system, 
it is convenient to normalize the variables, because in 
the vicinity of an hypothetical initial singularity physical 
variables would typically diverge, whereas at late times 
they commonly tend to zero 3]. Due to physical con¬ 
siderations normalization with the Hubble factor is an 
appropriate choice in cosmology. Besides, all available 
mathematical evidence suggests that Hubble-normalized 
variables are bounded into the past (that is, as the initial 
singularity is approached), and if there is a cosmologi¬ 
cal constant (or something that mimics it) it seems those 
variables will also be bounded into the future. Thus, even 
though the Hubble-normalized state space is unbounded, 
it is sensible to expect that the evolution equations will 
admit a past attractor and a future attractor. 


A. High energy limit 

In view of the discussion above, and following the usual 
practice, we use normalized (dimensionless) variables de¬ 
fined by 


0 

V6H' 

( 2 . 6 ) 

w 

(2.7) 

V ~ VSH' 

a 1 ~ n 

3 H 2 ' 

( 2 . 8 ) 

w= ^ . 

VS H 

(2.9) 


This coordinates will allow us analyzing the solutions 
of (l2. J2.2ll . and the cosmological models associated with 
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TABLE I: Location and existence conditions of the critical points of the HELA dynamical system 


Name 

X 

V 


W 


Existence 

W± 

0 

0 


±1 


All 7 and 

A 

x± 

±1 

0 


0 


All 7 and 

A 

XY±W+ 

7 

2A 


-7)7 

\/‘- 

7 

2A 2 

1 > A 2 > 

7 

2 

XY ± W- 

7 

2A 


-7)7 


7 

2A 2 

1 > A 2 > 

7 

2 

XY± 

A 

±vi- 

A 2 

0 


A 2 < 1 


0 

0 

0 


0 


All 7 and 

A 


them. As will become clear below, this choice of coordi¬ 
nates leads to a dynamical system which is best adapted 
for studying high energy limit solutions, so we will refer 
to it as our high energy limit adapted (HELA) dynamical 
system as opposed to a low energy limit adapted (LELA) 
that will emerge below after an alternative choice of co¬ 
ordinates. 

In addition, the variables will be related among them 
through 

S(w, x, y) + S n (w, x, y)z 1 ~ n = 1, (2-10) 

where S(w, x , y) = w 2 + x 2 + y 2 . 

The constraint imuii lets us “forget” about the evo¬ 
lution of one of the coordinates. Here we will choose 
the discarded coordinate to be 2 . Using the variables 
EH , equation Ell . and the conservation equations 
(IQ) and (IQ) we get the equations 

X ' = ^s ^ ^ n - ^ x2 ~ 2y2 + n ~ 2 ) w2 ) “ 


((n — 1) x ( 2x 2 + 'yxw 2 ) — 2A y 2 ) S) , 

(2.11) 

y'=^ (n( 2 x 2 +jw 2 ) - 


(2Ax + (n — 1) (2x 2 + jw 2 )) S ) , 

(2.12) 

u/=— (x 2 (2 n — 7) + 7 (n — 1 ) w 2 — 


7 y 2 — (n — 1) (2a; 2 + yw 2 )S) , 

(2.13) 

1 


If the real part of the three characteristic exponents is 
negative, the fixed point is asymptotically stable, i.e., an 
attractor. On the other hand, it is enough to have (at 
least) one characteristic exponent with positive real part 


where primes denote differentiation with respect to r e 
In a 3 . Equations 12.1112.1311 form the HELA dynamical 
system, which is defined on the phase space 

^ = {(x,y,w) : S < 1}. (2.14) 

To begin with the study of the evolution of our the 
dynamical system we have to find its fixed (or critical) 
points. The fixed points (x* ,y* ,w*) are given by the 
conditions 


x' (.x*,y*,w *) = 0 , 

(2.15) 

y'(x*,y*,w*)= 0 , 

(2.16) 

w' (a;*, y*, w*) = 0 . 

(2.17) 


In order to analyze the stability of the fixed points 
(x*,y*,w*) one studies the linearized dynamical sys¬ 
tem obtained by expanding equations TO EUt about 
those fixed points (see, e.g [j3])- Then one tries solutions 
in the form (x,y,w) = (ci, C 2 , C 3 ) e xt in the linear ap¬ 
proximation, and finds that the characteristic exponent 
A and the constant vector (ci, 02 , 03 ) must be respectively 
an eigenvalue and an eigenvector of the matrix 


/ dx' 

dx' 

dx' \ 

dx 

dy 

dw 


dtf_ 

dj_ 

dx 

dy 

dw 

duo' 

dw' 

dw’ 

\ Ite 

dy 

dw / 


to make the fixed point asymptotically unstable, i.e., a 
repeller. This repeller is a saddle point if at least one 
of the other characteristic exponents has a negative real 
part, in which case there is, apart from the unstable man- 
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ifold, a stable manifold containing the exceptional orbits 
that converge to the fixed point. The case in which the 
largest real part is precisely zero must be analyzed using 
other methods, for the linear analysis in unconclusive. 
The geometric form of the orbits near the fixed point is 
determined by the imaginary part of the characteristic 
exponents. If the three are real the fixed point is a node. 


A couple of complex conjugate exponents leads, except in 
degenerate cases, to an spiral point (the orbits are helices 
near the fixed point). In addition, when one of the ex¬ 
ponents is null the point is not hyperbolic and therefore 
structural stability cannot be guaranteed (the geometric 
form of the trajectories may change under small pertur¬ 
bations). 


TABLE II: Eigenvalues and dynamical character of the critical points of the HELA dynamical system 
(assuming n<l, 0<y<2 and A 2 < 1 and (3 = -y/(2 — 7 )( 4 y 2 + A 2 (2 — 9y)) ^ 


Name 


Eigenvalues 


Dynamical Character 


W± 


X± 


XY±W+ 


XY±W- 


XY± 



Q(7-2), 7(1 -n)) 

(2(1 -n), 1-J, IT A) 

(7 ~ 2)A - /3 (7 — 2)A +J3 \ 
h 4A 2 ’ 4A J 


non-hyperbolic for 7 = 0 , 2 , repeller for 0 < 7 < 2 (saddle for 0 < 7 < 2 ) 

non-hyperbolic for 7 = 2 or A = 1, always repeller 

non-hyperbolic for 7 = 0 , 2 , repeller for 7 > 0 
(saddle or spiral-out saddle for 0 < 7 < 2 ) 


(*->• 


non-hyperbolic for 7 = 0 , 2 , repeller for 7 > 0 
(saddle or spiral-out saddle for 0 < 7 < 2 ) 


(A 2 - 1, A 2 - 1, 2(1-n)A 2 ) 


non-hyperbolic for A 2 = 1 , 7 / 2 , always repeller 
(saddle for A 2 < 1 or A 2 < 7/2) 


O 


( 0 , 0 , 0 ) 


undefined 


Depending on the values of A and 7 the HELA system 
above has up to eleven critical points. We will denote 
them as W±, X±, XY±W_, XY±W+, XY±,0 and the 
first ten are closely related those found in Ref. [l0| , where 
the n = 0 case of our system was studied. 

Observe that the equations are invariant under the 
variable changes y —> —y, z —> —z, but not under 
x —> — x, and so in our numerical examples we will con¬ 
centrate on the region {x 2 + y 2 + z 2 > 1 , — 1 < x < 1 , y > 
0, w > 0}, which is a fourth of the unit sphere. This is 
equivalent to saying we are just considering expanding 
universes (H > 0). However, an analytic description of 
all eleven points is presented in the lines below and in ta¬ 
bles © and O- We will also set restrictions 0 < 7 < 2 
and A < 1 so that neither the barotropic fluid nor the 
scalar field have supraluminical sound speeds and the 
fluid satisfies the weak energy condition. We will also 


assume n < 1, as this is the case of interest. 

Although the location of the critical points of this dy¬ 
namical system does not depend on n .23|, the same is 
not true for their dynamical character. Note as well that 
all the critical points but the point O are located on the 
hypersurface 5 = 1, that is, they correspond to z = 0 
and H = 00, which in turn implies those are solutions 
to be obtained under the high energy limit approxima¬ 
tion. Interestingly, the fact that H = 00 means those are 
singular (big-bang) cosmologies. 

In what follows, and in order to complete the infor¬ 
mation provided in the tables, we will characterize the 
cosmological models represented by the fixed point liv¬ 
ing in the above mentioned fourth of the unit sphere. 

The first point, called W+, represents a solution com¬ 
pletely dominated by the fluid. The unstable character 
of these solutions agrees with what one might have antic- 
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ipated, are they are only expected to be relevant at early 
times. 

The second point, called X + , represents a solution 
completely dominated by the scalar field, more specifi¬ 
cally by its kinetic energy. Figs. n anc 10 s t ron gly suggest 
it is a past attractor. 

The third point, called XY+W+, depicts a tracking 
solution, neither the fluid nor the scalar field dominate 
completely 0. 

The fifth point, called XY + , represents a scalar field 
dominated solution, which is inflationary if A 2 < 1/3 
:20j. Figs. 0 and 0] show it is a unstable saddle in the 
asymptotic future. 

As for the sixth point, called O the linear analysis is 
not sufficient for determining its dynamical character, be¬ 
cause it has three null eigenvalues. One possibility is 
resorting to numerical inspection. As shown in the fig¬ 
ures below, this approach tell us the fifth point is the 
late time attractor of the system. Another tricky aspect 
of the analysis of this point is the sort of solutions it 
represents. Clearly, in the vicinity of (0,0,0) we have 
solutions in which cj) <C H, V <C H and p 7 <C H 2 , but 
nothing can be said about the ratios </> 2 /V , U/p 7 and 
(p 2 /p-y. This means we cannot say whether the solution 
represented by O is scalar field dominated, fluid domi¬ 
nated or tracking. Note as well that in the vicinity of 
O we have S n (w,x,y)z 1 ~ n — 1, which in turn implies 
3 H 2 — crp™ ot , that is, the points in that region represent 
solutions in which the corrections terms in the Friedmann 
equation dominate, and this obviously happens in the 
low energy regime. We have seen this system is good for 
the description of the points far away from that region 
(i.e. associated with solutions in the high energy regime), 
so we say the system is adapted to a high energy limit 
description and thus dub it HELA system, as opposed 
to the system adapted to a low energy limit description 
which will be used in the next subsection, which we dub 
LELA system. By reformulating the whole analysis using 
the LELA system as an alternative to the HELA one, we 
will be able to unveil the fine structure of the O point, 
and more importantly we will be allowed to tell what sort 
of solution is the late-time attractor. 



FIG. 1: Phase-space trajectories lying on the y = 0 hypersur¬ 
face for 7 = 1, n = 0.2 and A = 1. 



FIG. 2: Phase-space trajectories lying on the w = 0 hyper¬ 
surface for 7 = 1, n = 0.2 and an arbitrary A. 



B. Low energy limit 

We consider now dimensionless variables defined by 


x = (V3H^j " T, 

(2.19) 

y = " VV, 

(2.20) 

w = cr 2^ (y/3 H^j n y/p, 

(2.21) 

z = a 2n ^~ n i ^V3, 

(2.22) 


FIG. 3: XY projection of some phase-space trajectories for 
7 = A = 1 and n = 0.2. All trajectories cross the hypersurface 
w = 0.69. 


S(w,x,y) n + S(w,x,y)z 2(1 n) = 1. (2.23) 

Using variables ( 12 . 1912. 2 211 , the relation ( 12 . 2311 . and the 
field equations we get the equations 

x' =—!-— (2 (n — 1) (S n — 1) x 3 — 2 nxy 2 + 

2 nS v 

w 2 x ((1 + (n — 1) S n ) 7 — 2 n) + 


where the above variables are related by 
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FIG. 4: XW projection of some phase-space trajectories for 
7 = A = 1 and n = 0.2. All trajectories cross the hypersurface 
w = 0.69. 


2ny/Sy/l - S n y 2 A) , 




(2.24) 

y 

2 nS 

((1 + (n - 1) S n ) (2 x 2 + 


7 w 2 

) — 2nXx\/ S Vl — 5™^ , 

(2.25) 

1 

2 nS 

((n - 1) {S n - 1) w 3 7 - 


w (n 

l { y 2 -l)-z 2 ( 2(1 + (n 

-1)5”)))) (.2.26) 


Equations 12.2412.2611 form a dynamical system defined 
on the state space 

V = {{x,y,w) :S(w,x,y) < 1 } . (2.27) 

As discussed above this will be our LELA system, and 
like in the HELA one, the location of the critical points 
of this dynamical system does depend on n, but their the 
dynamical character does. 


The LELA system has always seven critical points 
W ±, X±, Y±, O, no matter the values of A and 7 . Once 
again, given the invariance under the transformations 


FIG. 5: Phase-space trajectories lying on the y = 0 hypersur¬ 
face for 7 = 1, n = 0.2 and A = 1. 



FIG. 6 : Phase-space trajectories lying on the w = 0 hyper¬ 
surface for 7 = 1, n = 0.2 and an arbitrary A. 



FIG. 7: XY projection of some phase-space trajectories for 
7 = A = 1 and n = 0.2. All trajectories cross the hypersurface 
w = 0 . 1 . 


TABLE III: Location and existence conditions of the critical 
points of the LELA system 


Name 

X 

y 

Z 

Existence 

x ± 

±1 

0 

0 

All 7 and A 

Y ± 

0 

±1 

0 

All 7 and A 

W ± 

0 

±1 

0 

All 7 and A 

O 

0 

0 

0 

All 7 and A 


y —> —y, z —> —z, combined with the lack of invariance 
under the transformation x —> —x, we will concentrate on 
the region {x 2 +y 2 + z 2 > 1 , — 1 < x < 1 , y > 0 , w > 0 } in 
our numerical examples, and as explained above this will 
mean we will address expanding universes only. However, 
an analytic description of all seven points is presented in 
the lines below and in tables ITTTI and ITVl We will also set 
the same restrictions on 7 , A and n as above. 

We have seen that the HELA system is not suitable for 
the description of solutions satisfying 3 H 2 « ap n . Un¬ 
der the definitions used to construct our LELA system 
such solutions live in the region S ss 1 , so the solutions 
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FIG. 8 : XW projection of some phase-space trajectories for 
7 = A = 1 and n = 0.2. All trajectories cross the hypersurface 
y = 0 . 01 . 

represented by the fixed points of W±, X± and Y± of 
our LELA system are also solutions where corrections 
are important. From lt2~23ll we see those points are char¬ 
acterized by 2 = oo, which is a condition equivalent to 
H = 0. Since, in an expanding universe H goes to zero 
at late times, the previous conclusion is consistent with 
the requirement that at late times acceleration appears 
due to the domination of corrections. 

In contrast the fixed point O represents solutions for 
which 3 H 2 ft p, that is, solutions in a regime in which 
corrections are negligible, but since these solutions were 
already investigated during our analysis of the HELA 
system [24|, we will not bother to study the features of 
this O point. 

In such limit, the critical points are X±, Y±, and 


W ±, In (IIIII) we summarize the information about these 
points. 

The first point, called X + , is located at (x,y,w) = 
(1,0,0). It is saddle if 7 7 ^ 2 and repeller otherwise. It 
represents a solution completely dominated by the scalar 
field, and more specifically by its kinetic energy. 

The second point, called Y+, is located at (x,y,w) = 
(0,1,0). Its eigenvalues are undefined because when the 
point is approached from different directions different re¬ 
sults are obtained. Numerical analysis suggests it is the 
late time attractor, it represents a solution completely 
dominated by the scalar field, and more specifically by 
its potential energy. 

The third point, called W + , is located at (x,y,w) = 
(0, 0,1). It is a saddle if 7 7 ^ 0, 2, a repeller if 7 = 2 and it 
is undefined if 7 = 0. It represents a solution completely 
dominated by the fluid. 

Finally the fourth point, called O is located at 
(x,y,w) = (0,0,1). Its dynamical character is also un¬ 
defined, but numerical analysis indicates it is a repeller. 
The choice of variables leading to the LELA system does 
not allow for an investigation of this fixed point, but it 
must be noted that it comprises all the points in the 
HELA system except the O point, so there is closes the 
circle, because the character of those points was investi¬ 
gated in the previous section. 

Note that the problems to define the dynamical char¬ 
acter of some points in the HELA and the LELA systems 
should not come as a surprise, because this is a common¬ 
place problem in points lying outside the variety (points 
at which the equations or their derivatives blow up). 


TABLE IV: Eigenvalues and dynamical character of the critical points of the LELA dynamical system 

(assuming n < 1, 0 < 7 < 2 and A 2 < 1) 


Name 

Eigenvalues 

Dynamical Character 

x± 

(- 2(1 -n), 1 , 1 -J) 

non-hyperbolic for 7 = 2 , always repeller (saddle for 0 < 7 < 2 ) 

Y± 

undefined 

undefined 

W± 

(_(l_ n ) 7 ,_l + l, 2 ) 

non-hyperbolic for 7 = 0 , 2 , repeller for 7 > 0 (saddle for 0 < 7 < 2 ) 

0 

( 0 , 0 , 0 ) 

undefined 


III. COSMOLOGICAL CONSEQUENCES standard cosmology, for the description of Cardassian 

models we find that the first complication stems from 


The cosmological consequences of this analysis are 
simple but important. As compared to the situation in 

















the necessity of introducing an additional variable, which 
we call z. 

Our numerical analysis tell us that the past attractors 
correspond to S(w,x,y ) = 1, and because of the con¬ 
straint the latter enforces z 1 ~ n = 0 which with in turn 
implies the recovery of the usual form of the Friedmann 
equation. In the case of models expanding from an initial 
singularity will and for (n < 1) we then conclude that the 
past attractors corresponds more specifically to to z = 0, 
and from the definition of z we see that those are so¬ 
lutions with an initial singularity. Summarizing, from 
the perspective of dynamical systems Cardassian mod¬ 
els with a fluid and a scalar field with an exponential 
potential will preferably have a big bang. 

More specifically, the early-time attractor is a solu¬ 
tion completely dominated by the kinetic energy of the 
scalar field and satisfying p oc a~ 6 , and its evolution is 
indistinguishable from that of perfect fluid models with 
a p = (7 — 1 )p equation of state under the identifica¬ 
tion n = 7/2, and the condition for inflation is simply 
n < 1/3. 

In contrast, the late-time solution attractor is a solu¬ 
tion completely dominated by the potential energy of the 
scalar field and satisfying p = const., that is, the de Sit¬ 
ter solution, and since H = const, too, inflation proceeds 
for any value of n. 

Interestingly, there are no tracking late-time attrac¬ 
tors, this is an important difference with respect to the 
behavior in standard (non-Cardassian) models. 

IV. CONCLUSIONS 

Cardassian models have been proposed as yet one more 
possible explanation for late-time acceleration. The main 
interest of the proposal is it involves only matter and ra¬ 
diation and does not invoke either vacuum energy or a 
cosmological constant. The idea consists in introducing 
a modification to the Friedmann equation, so that the 


effects of the modification become important at low red- 
shift. 

We have concentrated here on modifications of the 
form 3 H 2 3 * — p oc p n with n < 1, and we have studied 
its asymptotic behavior assuming p is made up of two 
contributions: the energy density of a perfect fluid with 
a p = (j — l)p equation of state and a self-interacting 
scalar field with an exponential potential. 

Our analysis falls mainly on the analytical side, but we 
have also carried out some numerical investigations. We 
constructed two dynamical system arising from the evo¬ 
lution equations. The first system is best suited to telling 
us about the asymptotic behavior in the high energy limit 
and we have called it our HELA system. The The sec¬ 
ond system, in contrast is best suited for a description of 
the asymptotic behavior in the low energy limit and we 
have called it our LELA system. We have found all the 
fixed points and we have characterized dynamically most 
of them by analytical means. However, the information 
regarding the dynamical character of some of the most 
important points seem to be only accessible by numerical 
methods. 

Our analysis allows us to say that the late-time attrac¬ 
tor of the system is the de Sitter solution, and whatever 
the value of n, there will only be inflation. We have also 
found that the early-time attractor is a solution with a 
big-bang, and that there are not tracking late-time at¬ 
tractor. Basically, the behavior in the high energy regime 
is the same as in standard cosmology but the modifica¬ 
tions make the structure of late-time attractors rather 
different. 
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